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ABSTRACT 

In an attempt to detect cosmic microwave background (CMB) anisotropy on arcmin 
scales, we have made an 8. 7- GHz image of a sky region with a resolution of two arcmin 
and high surface brightness sensitivity using the Australia Telescope Compact Array 
(ATCA) in an ultra-compact configuration. The foreground discrete-source confusion 
was estimated from observations with higher resolution at the same frequency and 
in a scaled array at a lower frequency. Following the subtraction of the foreground 
confusion, the field shows no features in excess of the instrument noise. This limits 
the CMB anisotropy flat-band power to Q/iat < 23.6 fiK with 95 per cent confidence; 
the ATCA filter function* Fi in multipole /-space peaks at l^ff = 4700 and has half 
maximum values at I = 3350 and 6050. 

Key words: cosmic microwave background - cosmology: observations - techniques: 
interferometric. 



1 INTRODUCTION 

Anisotropies in the cosmic microwave background (CMB) 
are usually described in terms of the coefficients Cim of 
their spherical harmonic decomposition; Ci = {| Cim \'^) 
represents the anisotropy power at multipole order I. Sachs- 
Wolfe anisotropies in a scale-invariant matter power spec- 
trum give flat l{l + l)Ci and any experiment sensitive to 
CMB anisotropies in an /-space window Fi may be expected 
to measure a sky temperature variance 

where Qfiat denotes the quadrupole normalization of a flat 
CMB anisotropy spectrum. The observing scheme and data 
analysis procedures determine the window function Fi for 
any experiment. The implications of any observed CMB 
temperature variance (AT)^ may be quoted in fairly model 
independent terms by the inferred Qfiat- This Qfiat may be 
viewed as being a measure of the CMB anisotropy power 
spectral density within the /-space window defined by Fi. 

Recent observations of the anisotropy in the CMB on 
large and intermediate angular scales (/ ^ 500) have begun 
to provide interesting constraints on theories of structure 
formation and the parameters of cosmological models (see, 
for example, Lineweaver and Barbosa (1997) for an updated 

*The ATCA filter function is available at the website 
www. atnf. csiro.au/research/ cmbr/cmbr^atca.html 



compilation of data points and a discussion of their impli- 
cations). 

Primary CMB anisotropies at small angular scales (/ ^ 
500) are expected to be relatively damped in most struc- 
ture formation models owing to the the thickness of the last 
scattering surface and the diffusion damping of sub-horizon 
scale baryon fluctuations in the pre-recombination epoch. 
However, flat-band powers comparable to the Q fiat detected 
by COBE DMR (Bennett et al. 1996) may be predicted at 
multipoles / ^ 500 in certain cosmological-constant (A) and 
open-universe manifestations of baryon isocurvature mod- 
els (Hu & Sugiyama 1994; Hu, Bunn & Sugiyama 1995). 
The anisotropy power at large / may be critically dependent 
on the reionization history: small-scale anisotropies may be 
suppressed by early reionization. 

Secondary anisotropies may be generated at the later 
last-scattering-surface in a reionized universe predominantly 
owing to second-order mode coupling between density per- 
turbations and bulk velocities (the 'Vishniac effect'); this 
could significantly contribute to arcmin-scale (/ J5 10'^) 
anisotropies particularly if the ionization fraction is high at 
late times (Hu, Scott & Silk 1994). Persi et al. (1995) esti- 
mate that this second-order Doppler effect may contribute a 
band power Q fiat ~ 0.2 fiKa^tl ^ 10^ in CDM and CDM-I-A 
universes and that the anisotropy power may be an order of 
magnitude higher in certain baryon isocurvature models. 

Jones et al. (1997) report the detection of a ~ 100 fiK 
negative feature in a 2-arcmin resolution image of a sky 
patch that has no obvious cluster of galaxies along the line 
of sight in either optical or ROSAT X-ray images. A sen- 
sitive low-resolution image of a 'blank' field with the VLA 
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(Richards et al. 1997) is also reported to show a negative fea- 
ture approximately 25" x 65" in size and with a peak central 
decrement of —250 /itK. It has been suggested that these may 
be CMB decrements owing to the inverse-Compton scatter- 
ing of CMB photons, the Sunyaev-Zeldovich (S-Z) effect, in 
distant concentrations of hot gas. It may be noted that Hat- 
tori et al. (1997) recently reported the discovery, in its X-ray 
emission, of a hot gas concentration at a redshift z w 1 whose 
properties arc similar to the gaseous halos in rich clusters 
of galaxies; however, the 'cluster' appears to have only one 
visible galaxy. 

Anisotropics arising from the S-Z effect, in a cosmo- 
logical population of groups and clusters of galaxies con- 
taining hot intra-clustcr gas may be an important cause of 
anisotropy power at I ^ 10"^ (Persi et al. 1995; Bond & My- 
ers 1996). The S-Z anisotropies are generically non-Gaussian 
and are expected to have phase-correlations between dif- 
ferent I modes; their statistical description requires higher- 
order correlations apart from the Ci power spectra. Hot gas 
in groups and clusters in scale-invariant CDM universes nor- 
malized to give CT8 = 1 (the magnitude of dynamical cluster- 
ing at the present time is quantified by as which is the rms 
mass fluctuations in 8/i~^-Mpc spheres) is expected to con- 
tribute band powers Q fiat ~ 2 /jK at I ^ 1-5 x 10^ (Bond & 
Myers 1996). In COBi^normahzed tilted CDM, MDM and 
models in which the shape factor T = Qoh is approximately 
0.2, the predicted ag agrees better with galaxy clustering 
and in these models the band powers at large I are 1-2 orders 
of magnitude smaller. The dominant contribution at high I 
may come from quasar-ionized hot gas bubbles (Aghanim et 
al. 1996), whose existence is based on plausible theoretical 
inferences. 

Computations of the expected secondary anisotropies 
involve non-linear gravitational dynamics and hydrody- 
namic simulations and are critically dependent on the ther- 
mal history of the gas that may in turn depend on the astro- 
physical evolution in populations that cause the ionization. 
It follows that observations of the arcmin-scale anisotropy, 
corresponding to measurements of the CMB anisotropy at 

1 ^ 10^, could constrain the structure formation theories. 

2 OBSERVATIONS WITH THE AUSTRALIA 
TELESCOPE 

The advantages of Fourier-Synthesis imaging telescopes, and 
in particular the design features of the Australia Telescope 
Compact Array (ATCA; see The Austraha Telescope 1992) 
that make it specifically advantageous for high-brightness- 
sensitivity imaging, were detailed in Subrahmanyan et al. 
(1993) and we restrict ourselves to giving a synopsis of the 
methodology here. Our observing strategy has been to make 
full Earth-rotation synthesis observations of 'empty' fields 
in a special ultra-compact 122-m array configuration - with 
five 22-m diameter antennae located 30.6 m apart in an E-W 
line - in order to maximize the brightness sensitivity of the 
imaging. The three baselines between antennae spaced 61 m 
apart are used to construct a model of the foreground con- 
fusion; this is then subtracted from all of the visibility data 
and the four baselines between antennae 30.6 m apart are 
used to synthesize an image with high brightness sensitiv- 
ity. Because the confusing sources are measured simultane- 
ously and using baselines between the very same antennae, 



errors in the estimation of confusion owing to variability in 
the foreground sources and calibration errors are eliminated. 
Deconvolution errors are avoided because the confusion is es- 
timated as a model fit to visibility data. These observations 
were made at the highest available frequency of 8.7 GHz to 
minimize discrete source confusion. 

The field was separately observed in a nearly scaled 244- 
m array at 4.7 GHz, with five ATCA antennae spaced 61-m 
apart along the E-W line, to examine the spectral indices of 
any features identified in the region as foreground sources 
and ensure that the sources subtracted indeed have spectra 
consistent with optically-thin thermal or synchrotron emis- 
sion. If the 4.7-GHz images are made with the same sensi- 
tivity to point-source flux density as the 8.7-GHz images, 
the lower-frequency images that are made with the same 
angular resolution would be less sensitive to CMB temper- 
ature fluctuations by a factor (4.7/8.7)", where a = 2 is 
the spectral index of the CMB anisotropy at these frequen- 
cies. However, the 4. 7- GHz images would have a relatively 
greater sensitivity to extended synchrotron emission by a 
factor (4.7/8.7)", where a ^ —0.7 is now the steep spectral 
index of extended synchrotron sources. To summarize, the 
4.7-GHz image could potentially reveal extended foreground 
sources that may be resolved, and therefore absent, in the 
61-m baseline data obtained at 8.7 GHz. 

Sky regions - selected using radio source catalogues to 
be relatively devoid of bright sources - were first imaged at 
20-cm wavelength with the ATCA. An examination of these 
images led to the selection of a field centre which had no 
sources with fiux densities exceeding 1 mjy (at 20-cm wave- 
length) within 7 arcmin radius (the ATCA antennae have 
the first null of the 8.7-GHz primary beam at 7 arcmin). 
The J2000.0-epoch coordinates of the field centre arc RA: 
03*^ 16"^ 26^:00, DEC: -49° 47' 57'.'00. The field has been 
chosen to be located close to declination —50° so that the 
projected antenna spacing is as close as possible without 
shadowing to the 22-m antenna diameter at large hour an- 
gles: this maximizes the brightness sensitivity of the array. 

The field was originally observed in 1991 July-August 
and December in the 122-m array and using a pair of 128- 
MHz bands centred at 8640 and 8768 MHz; results of the ob- 
servations made during these periods were reported in Sub- 
rahmanyan et al. (1993). Subsequently, the 8.7-GHz front- 
end amplifiers were changed to HEMT devices and the sys- 
tem temperature at this frequency improved from about 75 
to 43 K. The field was reobserved in 1994 July in the 122- 
m array using the same pair of bands. The total effective 
observing time obtained in this array now corresponds to 
about 50 h integration with a system temperature of 43 K 
and two 128-MHz bands. The rms sensitivity of the imaging 
has improved by a factor 1.6 as a result of the additional 
observations made with the improved receivers. All obser- 
vations were made in dual polarizations. The array phase 
centres during the observations were ofi'set about 1° from 
the antenna pointing centres so that imaging artefacts that 
often appear at the phase centre may be well removed from 
the sky region of interest. Off-line, the visibility data were 
phase-corrected to align the array phase centres with the 
field centre positions, calibrated in amplitude, phase and for 
the band-pass response, averaged in frequency over a useful 
band of 112 MHz. The flux-density scale was set by adopt- 
ing values of 2.84 and 2.79 Jy respectively for the primary 
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calibrator source PKS 1934—638 at the frequencies 8640 and 
8768 MHz (Reynolds 1994). 

In 1995 March, the field was observed in a 244-m ar- 
ray at 4800 MHz. The observing and calibration procedure 
adopted was the same as for the 122-m array observations; 
the flux-density scale was set by adopting a value of 6.22 Jy 
for the primary calibrator PKS 1934—638 at this frequency. 

All images shown below were made using visibilities 
with 'natural' weighting so as to obtain a high signal-to- 
noise ratio, the gridding, Fourier transformations and de- 
convolution, if attempted, were done using the aips routine 
IMAGR. All displayed images are centred at the coordinates 
of the field centre and, unless explicitely stated, have not 
been deconvolved. All images shown are in Stokes I and have 
not been corrected for the attenuation owing to the primary 
beam. The locus of the first null in the primary beam pat- 
tern at 8704 MHz, located at a radius of 7 arcmin, is shown 
in all the images as a dot-dashed circle. 

3 ANALYSIS OF THE OBSERVED FIELD 

An image of the field made at 4.8 GHz using the 61 and 
122-m baselines of the 244-m array observations is shown in 
Fig. 1. The image has been deconvolved and has a resolu- 
tion of 1.4 arcmin. Within a radius of 3 arcmin, correspond- 
ing to the half-power radius of the ATCA primary beam at 
8.7 GHz, three prominent sources are apparent and all these 
are detected with peak flux density exceeding 10 times the 
image thermal noise. 

In Fig. 2 we show an image of this field at 8.7 GHz that 
has been constructed using just the four 30.6-m baselines. 
This image has not been deconvolved, but clearly shows 
three peaks at the positions of the three sources apparent 
in the 4.8-GHz image. It may be noted that the synthesized 
beam sidelobes are very large for this image because of the 
poor u, V coverage (ordy a single spacing - 30.6 m - has been 
used in the imaging). 

The low-resolution 8.7-GHz image of the field (Fig. 2) 
has a high surface brightness sensitivity, but is confusion 
limited owing to discrete sources in the field. We estimate 
the foreground confusion owing to discrete sources in this 
field using the data obtained in the longer (> 30.6 m) base- 
lines. However, we have restricted ourselves to using just 
the 61-m baseline (omitting the 92 and 122 m baselines) for 
deriving the confusion model so that the synthesized beam 
is only a factor of two smaller as compared to that for the 
high-surface-brightness image (Fig. 2) from which the co- 
fusion is to be subtracted. This may ensure that confusion 
structures on scales up to about an arcmin, if present, may 
be included in the confusion model. 

We next show a deconvolved image of the field made 
at 8.7 GHz using the 61-m baselines (Fig. 3). Three sources 
appear once again, we estimate their positions and flux den- 
sities at 8.7 GHz from this deconvolved image and use these 
parameters as an input model while fitting a 3-component 
model to the 61-m visibility data. The fit estimated the fiux 
densities of the components to be 274, 182 and 110 /uJy. The 
primary beam at 8.7 GHz is expected to have attenuated the 
source intensities by factors 0.57, 0.39 and 0.51 respectively. 
Correcting for this attenuation, the components ought to 
have flux densities 481, 467 and 216 /nJy respectively. We 
have fit the 4.8-GHz visibilities measured using the 122-m 
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Figure 1. Deconvolved image of the field at 4800 MHz made 
witfi a beam of 96 X 72 arcsec'^ at a position angle (p. a.) of 52°. 
Contours at 27 fijy beam^i X (-3, 3, 4, 6, 8, 12, 16, 24, 32). 
The half-maximum size of the synthesized beam is shown in the 
bottom right corner as a filled ellipse. In this figure, as also the 
following three figures, the locus of the first null in the primary 
beam pattern at 8704 MHz is shown as a dot-dashed circle. 
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Figure 2. An image of the field at 8704 MHz made using just the 
30.6-m baselines. The resolution is about 2.2 arcmin. Contours at 
22 /^Jy beam-l X (-10, -8, -6, -4, 4, 6, 8, 10, 12, 14, 16). 

baselines (equivalent in angular resolution to the 61 m base- 
lines at 8.7 GHz) to a 3-component model that has compo- 
nent positions fixed at the locations of the sources in the 
8.7-GHz, 61-m visibilities: the flux densities of these compo- 
nents are estimated to be 790, 335 and 262 //Jy respectively. 
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Figure 3. A deconvolved image of tlic field at 8704 MHz made 
using just tfie 61-m baselines. The image lias a beam FWHM of 
71 X 57 aresee^ at a position angle (p. a.) of 0°. Contours at 25 
^Jy beam^i X (-4, -3, -2, 2, 3, 4, 6, 8, 10). The half-maximum 
size of the synthesized beam is shown in the bottom right corner 
as a filled ellipse. 
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Figure 4. The residual image of the field at 8704 MHz follow- 
ing the subtraction of the confusion model; this image has been 
made using just the 30.6-m baselines. The resolution is about 2.2 
arcmm. Contours at 22 /ijy beam"! X (-3, -2, -1, 1, 2, 3). 



do not have spectral indices a = 2 that would be expected 
of CMB anisotropics at these frequencies. 

The three-component model for the foreground confu- 
sion, derived from the fit to the 60-m visibilities obtained at 
8.7 GHz, was subtracted from all of the 8.7-GHz visibility 
data. Following subtraction of the confusion, wc constructed 
an image of the field at 8.7 GHz using just the 30.6-m base- 
lines (Fig. 4). No residual features are apparent in this im- 
age. Wo weighted the image pixel intensities by the primary 
beam pattern and determined the weighted rms in the pri- 
mary beam region of the field to be 21.2 /jjy beam~^. The 
image pixel variance appears consistent with that expected 
from the telescope thermal noise (see section 4.1), and in 
sections 5 and 6 below wo use this residual image to derive 
limits on random-phase CMB anisotropy on the 2-arcmin 
scale corresponding to the resolution in this image. 

4 CONTRIBUTORS TO THE IMAGE 
VARIANCE 

The processes which lead to the synthesized image are 
shown schematically in Fig. 5. The telescope thermal noise, 
foreground discrete sources and possibly CMB anisotropy 
contribute to the variance in the synthesized image. CMB 
anisotropy and foreground radio sources on the sky are 
attenuated by the primary-beam pattern of the antennae 
and convolved by the synthesized beam; their contribution 
may equivalently be viewed as a visibility-domain spatial- 
frequency filtering of the sky signals as viewed through the 
primary beam. The telescope thermal noise may be consid- 
ered to be additive Gaussian random fluctuations, that are 
characterized by an appropriate angular coherence function, 
but whose distribution on the sky image has no dependence 
on the primary beam. 

The primary beam plays an important role in the signal 
path from sky sources to visibility data (or equivalently the 
synthesized image) and the /-space filter function that deter- 
mines the coupling of CMB power to image variance depends 
on the antenna aperture illumination function. Therefore, 
instead of modelling the primary beam in the sky plane, 
we have approximated the aperture illumination by a func- 
tion whose parameters characterize the central blockage and 
edge taper (James 1987) and have used parameter values 
that give an antenna power pattern that fits the measured 
pattern out to the second null. 

The number of independent pixels in the sky image 
made with just the 30.6-m baseline data, or equivalently the 
number of degrees of freedom in the variance estimate, may 
be approximately the number of independent visibility mea- 
surements along the 170 m long visibility track. The first null 
of the ATCA primary beam power pattern is at 7 arcmin ra- 
dius at 8.7 GHz, this implies that the aperture illumination 
has an autocorrelation function with an effective diameter of 
about 17 m and hence the number of independent visibility 
measurements is about 10. 



Correcting for the attenuation owing to the primary beam at 
4.8 GHz, we estimate the true flux densities of these compo- 
nents to be 934, 439 and 319 /iJy respectively. The derived 
spectral indices of the components are —1.1, 0.1 and —0.7 
(we define the spectral index a as Si, ~ i/™); the components 



4.1 Telescope system noise 

The image variance resulting from the telescope thermal 
noise has boon measured by two methods. First, by com- 
puting the rms over large sky regions well outside the pri- 
mary beam area where the only contribution to the image 
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Figure 5. The signal flow path. It may be noted that the point spread function (PSF) corresponding to the sky sources - the synthesized 
beam - may differ from the PSF associated with the telescope thermal noise. 



variance is expected to be the telescope thermal noise. Sec- 
ond, by separately imaging the visibilities obtained as the 
XX and YY correlations between the two orthogonal linear 
antenna signals X and Y and computing the rms over the 
(XX— YY)/2 image: the Stokes I component of sources in 
the field cancel in this combined image and if we assume 
that the Stokes Q, U and V flux densities are negligible 
compared to the thermal noise, the system thermal noise 
will be the only significant contributor to the variance in 
this combined image. The two methods have yielded consis- 
tent values and the rms thermal noise in the sky image is 
estimated to be 22.8 /xJy beam~^ in the image of the field 
made using the 30.6-m baseline data. Adopting an antenna 
efficiency of 0.64 at 8.7 GHz (James 1985) the rms noise in 
the image is consistent with the system temperatures be- 
ing effectively 72 and 43.5 K respectively during the 1991 
and 1994 epoch observations; the later value is in agreement 
with 1994 epoch measurements of the system temperature 
made with microwave absorbers placed over the feed horns 
(Gough 1994). 

Along the visibility tracks in the u, v-domain, the signal- 
to-noise ratio of the visibility data may vary because, for 
example, the total system temperature may systematically 
vary with antenna elevation and hence with position angle in 
the u, v-plane, or because data points may have difi^erent av- 
eraging times. Visibility data are almost never weighted pro- 
portional to the rms thermal noise; they are often weighted 
inversely with the thermal noise variance. Therefore, the im- 
age thermal noise is described by a coherence function that 
is usually not derivable from the synthesized beam. To com- 
pute the angular coherence function of the thermal noise, or 
equivalently the distribution of thermal noise power in the 
u, v-domain, we use images made well outside the primary 
beam area: we compute the 2D autocorrelation function of 



these sky images and perform a 2D Fourier transform to the 
u, w-domain. This provides us with a u, v-domain filter ap- 
propriate for the telescope thermal noise. Our simulations 

of the thermal noise component in the sky images have 
assumed the noise to be Gaussian random, with a power 
spectrum identical to this 2D u, u-plane filter function and 
normalized to give an image rms noise same as the above- 
mentioned measured values. 

The variance in the sky images have been determined 
by weighting the pixel intensities with the primary beam 
attenuation. We have simulated images with solely the ex- 
pected thermal noise contribution, computed the weighted 
image variances and obtained their distribution function. 
It may bo noted that because the variance estimate is the 
weighted sum of squares of random variables, the variance 
is not expected to follow a distribution. 



4.2 Foreground radio source confusion 

The radio sources in the sky are assumed to be Poisson dis- 
tributed and the differential source count N{S) is expected 
to be 



N{S) = 20.95" 



arcmin ^ /xJy 



(2) 



at 8.7 GHz, where S is the source flux density in jj,Jy 
(Partridge et al. 1997 and references therein). The ATCA 
Fourier-synthesis images have been constructed without any 
zero-spacing visibility components and at every image pixel 
position the expectation for the flux density is zero despite 
the presence of discrete sources in the sky. The expectation 
of the variance contribution from these sources at any image 
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position {60, (j>o) is given by 

y=yinax 

= J \b{e,(l>)\^d9d4> j y^N{y)dy, (3) 

y=ymiTi 

where y-min and ymax are the lower and upper hmits to the 
flux densities of sources present in the sky region and 0) 
represents a beam pattern that is the product of the pri- 
mary beam centred at the antenna pointing centre and the 
synthesized beam centred at {6o,4>o)- 

We reckon, based on our 1.4-GHz survey of the region, 
that the field has no sources exceeding about 0.5 mjy at 
8.7 GHz. The variance in the image made using the 30.6-m 
baseline data (Fig. 2) is a factor of two higher than that ex- 
pected owing to these sources assuming the Partridge et al. 
source counts. However, there is a reasonable (11 per cent) 
statistical probability that the 'excess' variance in the field 
is foreground confusion. The number of sources (3) detected 
with flux density exceeding 100 /LtJy - that constitute the 
model for confusion in the fleld - and their combined flux 
density, are also a factor of two larger than expected. A pos- 
sible explanation for these is that our 3-component model 
for confusion includes weaker sources that blend and give 
responses above the detection limit. 

Following the subtraction of the confusion model, the 
residual image shown in Fig. 4 does not show evidence for 
any significant 'excess' variance that may be owing to resid- 
ual weak foreground sources that have not been included 
in the model. For the assumed source counts, the expected 
variance contribution from sources weaker than 100 juJy is 
(50 /iJy beam"^)'^. The residual image shown in Fig. 4 has 
a variance of (21.2 /ijy beam^^)^ and the probability of ob- 
serving a sample image with a variance as low as this is 
2 X 10"* if sources below 100 /iJy have remained unsub- 
tracted! If we conjecture that the model derived from the 
61-m baseline visibility data has incorporated all sources 
above 25 /iJy, corresponding to the rms thermal noise in the 
image made with these data, we may then expect that the 
residual sources in the field that are weaker than this limit 
may contribute a variance of about (33 /iJy beam"^)'^ in im- 
ages made from the residual 30.6-m baseline data. Even in 
this optimistic scenario, the probability of observing a sam- 
ple with a residual variance as low as that in the observed 
field is as small as 0.3 per cent! 

To summarize, the image variance appears to be in ex- 
cess of the expectation prior to any subtraction of sources, 
but after the confusion model is subtracted the image vari- 
ance is below expectations that are derived assuming a 
minimal residual and unsubtracted confusion. We are led 
to believe that either the derived confusion model includes 
sources weaker than 25 /jJy, or perhaps that the faint sources 
are clustered. We are expanding our observations to several 
sky regions and expect to have a better understanding of 
the contribution from confusion following an examination of 
the residual fluctuations in a larger sky area. While deriving 
limits on CMB anisotropy in sections 5 and 6 based on our 
observations of the fleld, we assume that the residual im- 
age has contributions from the telescope thermal noise and 
CMB anisotropy alone. 

We have examined the possibility of using the skew co- 
efficient 03 as a possible indicator of residual foreground 



confusion. 03 is deflned as 

«3 = (■^^)3' ^itll "^3 = X]^?' "^2 = (4) 

and where the summations are over pixel intensities that 
are weighted by the primary-beam attenuations at the pixel 
locations. We find that the image in Fig. 2 has as = 0.13 
and the residual image in Fig. 4 has 03 = —0.03. Our simu- 
lations of foreground confusion show that the derived skew 
coefficient may be expected to have a standard deviation of 
as much as 0.57 about the mean skew, and that the expecta- 
tion for the skew coefficient is just 0.11 and 0.05 respectively 
if foreground sources below 100 and 25 /xJy remain unsub- 
tracted. Because of the limited number of independent sky 
pixels in the image of the field, the skew coefficient is not 
a useful statistic for characterizing the residual confusion; 
it may prove useful in statistical analyses of future observa- 
tions that cover multiple flelds. 

4.3 CMB anisotropy 

We model the CMB anisotropy as a Gaussian random fluc- 
tuation in the sky temperature that is completely described 
by its power spectrum coefficients Cj. The expectation for 
the image variance will then depend on the telescope filter 
function (TFF) Fi and the nature of the anisotropy power 
spectrum over the range in Z-space that corresponds to the 
pass-band of the filter. 

We assume that the anisotropy power spectrum is 'flat' 
in /(/ -|- l)-space and that the Cj coefficients are given by 

Ci = {Qnatf (^^(7^) ' (5) 

where Qfut denotes the quadrupole normalization of the 
power spectrum. The temperature variance in the sky image 

is then expected to be 

I 

and this leads to the expression in equation (1). 

We compute the image variance as a weighted mean of 
the squares of the pixel intensities (image pixel intensities 
are in units of Jy beam^^) using weights for the intensi- 
ties that are the primary-beam attenuations at the pixel 
positions. As discussed in Appendix A, we have computed 
the ATCA TFFs at the different image pixel locations us- 
ing equation (A4) and averaged them, using weights that 
are again the primary-beam attenuations, to obtain the fil- 
ter function Fi, in units of (Jy beam~^ K"^)^, appropriate 
to the computed image variance. We have computed the 
Fi's using the ATCA synthesized beams corresponding sep- 
arately to the images made with the 30.6 and 61-m base- 
line data, they are shown in Fig. 6. In the large-Z regime 
of the ATCA Fi , the expected image variance (in units of 
(Jy beam~^)2) is approximately given by 

^S' = {Qnatf Y^f^F. (7) 
I 

The TFF corresponding to the 30.6-m baseline data, 
shown as the continuous line in Fig. 6, is the filter function 
corresponding to the residual image shown in Fig. 4. It has 
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Figure 6. ATCA telescope filter function. The continuous line shows the TFF corresponding to the image made with the 30.6-m baseline 
data, the dashed line shows the filter corresponding to the 61-m baseline data. 



a peak of 0.48 at I = 4700 and has half-maximum values at 
I = 3350 and 6050. Summation over Z-space yields 

If we adopt a 'flat' CMB anisotropy spectrum with nor- 
malization Qfiat = 18 /iK corresponding to the detected 
power at multipoles I <, 20 (4-yr COBE DMR results in 
Bennett et al. 1996), the expected variance contribution in 
the 30.6-m baseline image is expected to be 18^ x 0.765 = 
(15.7 /xJy heam-^f. 

A Fourier synthesis telescope with finite-aperture ele- 
ments measures visibilities that may be considered to be the 
convolution of the all-sky visibility with the autocorrelation 
of the antenna aperture illumination. The visibility measure- 
ments may be viewed as samples along ID (u, i;)-tracks of 
this 2D (u, v)-plane visibility function. Because the ATCA 
antenna aperture extends 11 m in radius, its autocorrelation 
function will extend to a radius of 22 m. It may be noted 
that this autocorrelation function has a diameter of 44 m 
which exceeds the spacing between the ID {u, w)-tracks cor- 
responding to the 30.6 and 61 m baselines. Therefore, the 
visibility data along the 30.6 and 61-m baseline visibility 
tracks are individually averages over regions of the all-sky 
visibility function and these regions mutually overlap. This 
results in a partial overlap between the TFFs correspond- 
ing to the 30.6 and 61-m visibility data. At any instant, 
there are three 61-m baselines as compared to four 30.6-m 
baselines; therefore, the 61-m baseline image is less sensitive 
to unresolved foreground sources as compared to the 30.6- 
m baseline data, but by only a factor 0.87. The peak of the 
TFF corresponding to the 61-m baseline data is 0.087, a fac- 



tor 0.18 of the peak of the 30.6-m baseline TFF. Summing 
over the 61-m baseline TFF, J2l2Fi/{5l) = 0.085: CMB 
anisotropy with a flat spectrum will be expected to con- 
tribute a variance in the 61-m baseline image that is a factor 
0.11 of the variance in the 30.6-m baseline image. Adopting 
a normalization Qfiat = 18 ^K, the 61-m baseline image is 
expected to have a variance contribution (5.2 /xjy beam~^)^ 
from CMB anisotropy. Because (a) the two TFFs have only 
a small overlap and are largely sampling CMB power at dif- 
ferent multipole ranges and (b) the 61-m baseline TFF has 
a greatly reduced sensitivity to CMB anisotropy, the con- 
fusion model that is derived from the 61-m baseline data is 
unlikely to signiflcantly reduce the CMB anisotropy power 
expected in the residual 30.6-m bsiseline image. 

We have simulated sky images by generating 'flat' spec- 
trum CMB power in the {u, w)-plane: pixels of size 5u, 5v 
that are at distance y/u^ +v'^ from the centre of the (m, v)- 
plane are given Gaussian random complex conjugate visibil- 
ities with variance 

{ATf = {Qfiatf Sudv. (9) 

The (m, i))-plane CMB anisotropy model is then filtered by 
the TFF corresponding to the 30.6-m baseline data and 

inverted to form simulated sky images. We have thereby 
obtained distributions of the sample variance from simu- 
lated images that have variance contribution from CMB 
anisotropy alone. It may be noted here that the radial dis- 
tribution of power in the (m, u)-plane in the case of our 'fiat' 
model CMB anisotropy differs from that for thermal noise 
and, therefore, the two contributions differ in the number of 
degrees of freedom in their contributions to image variance. 
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Figure 7. Likelihood functions. The dashed Hne corresponds to the null hypothesis Hi, the continuous curve corresponds to the 
hypothesis IH2 that Qfiat = 22.4 ^K. The dotted line corresponds to the distribution function assuming that 'flat' spectrum CMB 
anisotropy with Q = 22.4 /iK is the only contributor to the image variance. 



5 LIMITS ON CMB FLAT-BAND POWER 

We use the likelihood-ratio test to derive limits on pos- 
sible CMB flat-band power in the data. We assume that 
the residual image has contributions only from the telescope 
thermal noise and possibly CMB anisotropy. The observed 
variance in the residual image shown in Fig. 4 is 



2 



(21.2 ^iJy beam"^)^ 



(10) 



We adopt the null hypothesis Hi that the variance con- 
tribution from CMB anisotropy is zero. 



Hi : acMB = 0. 



(11) 



Simulations of sky images that have only thermal noise with 
an expected variance (22.8 /xJy beam~^)^ then yields the 
probability distribution for the sample variance. We show 
this distribution as the dashed line in Fig. 7 and represents 
the likelihood function P(a^jj^|0) which is the probability of 
obtaining any observed variance conditional on Hi. 

We next hypothesize that the sky has flat-band CMB 
anisotropy quantifled by the normalization Qfiat, «-e., 



H2 : o-cMB corresponding to Q 



flat- 



(12) 



We simulate sky images with purely CMB anisotropy to ob- 
tain the distribution function for its variance contribution; 
wc show this distribution in Fig. 7 as a dotted lino for the 
specific choice Qfiat = 22 A /xK. The distribution function 
for the observed image variance, with contributions from 



thermal noise and CMB anisotropy, will be a convolution of 
the two individual probability distributions. This distribu- 
tion, shown as a continuous curve in Fig. 7, represents the 
likelihood function V{(Joba\Q fiat) which is the probability of 
obtaining any observed variance conditional on H2 . 

The size a of the test is the probability of rejecting 
H2 when it is true (a type I error). For the choice Qfiat = 
22 A fiK, a = 0.05 and, therefore, Qfiat = 22 A //K is a 95 
per cent confidence upper limit. 

The power /3 of the test is defined as: /3 = 1 — proba- 
bility of accepting H2 when Hi is true (a type 11 error). For 
our choice of Qfiat = 22.4 /iK, the power /3 = 0.47. In order 
to increase the power of the test to a value /3 = 0.5, we may 
increase the decision variance to a^f,^ = (21.8 /xJy beam~^)^ 
and change the hypothesis H2 to correspond to a choice 
Qfiat ~ 23.6 /iK. For this choice of parameters, the tost will 
reject H2 with size a = 0.05 (95 per cent confidence) and 
the test will simultaneously have a power /3 = 0.5. 

To summarize, the ATCA observations of the field place 
an upper limit of Qfiat < 23.6 /uK with 95 per cent confi- 
dence in an /-space filter that peaks at I — 4700. The filter 
has half-maximum values at I = 3350 and 6050. 

6 LIMITS ON CMB ANISOTROPY WITH A 
GAUSSIAN-FORM AUTOCORRELATION 
FUNCTION 

The CMB anisotropy is sometimes assumed to have a 
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Table 1. ATCA limits on GACF anisotropy models. 

Coherence scale Size of the test Power of the test CMB temperature rms AT/T 

1 i 





a 


/3 




CS /To 


0.'5 


0.05 


0.44 


58 fj,K 


2.1 X 10-5 


0.'5 


0.05 


0.50 


63 ^iK 


2.3 X 10-5 


l.'O 


0.05 


0.44 


42 /iK 


1.5 X 10-5 


l.'O 


0.05 


0.50 


45 ^iK 


1.6 X 10-5 


2.'0 


0.05 


0.44 


63 AtK 


2.3 X 10-5 


2.'0 


0.05 


0.50 


68 ^lK 


2.5 X 10-5 



To = 2.73 K 



Gaussian autocorrelation function (GACF). The sky tem- 
perature is modelled to have a Gaussian distribution with 
zero mean and the autocorrelation of the sky temperature 
is taken to be of the form 



C(0 = Co 



(13) 



Parameter Co ~ {T^{9,(f>)) represents the variance of the 
sky temperature and parameter ~ (— C(0)/C"(0))^ rep- 
resents the coherence scale (C" (5) denotes the second deriva- 
tive of the autocorrelation function). Assuming that the co- 
herence scale is a small angle, the power spectrum of the 
CMB temperature anisotropy may be approximately repre- 
sented as i-space coefficients 



Cl = 2'JtCo^c 



(14) 



This model may be expected to contribute an image variance 
(ASf = iCod ^(2/ + l)e"(~) Fi 



(15) 



to the ATCA residual sky image. In this expression, is 
the coherence scale in radians and Fi is the ATCA TFF 
corresponding to the variance estimate in the residual image 
of the observed field. 

Assuming that the sky temperature variance Co is an 
invariant, the expected image variance (AS)^ is a maximum 
for an anisotropy model with = 1 arcmin. We have sim- 
ulated sky images with the CMB anisotropy modelled to 
have a GACF and derived the distribution function for the 
image variance contribution. Limits to the model parameter 
Co may be derived for specific choices of the coherence scale 
5c using the likelihood ratio test described in the previous 
section. 

Adopting a value ~ 1 arcmin, we find that the ob- 
served image variance of a^^,^ — (21.2 /ijy beam-^)^ implies 

1 

a 95 per cent confidence upper limit of 42 /xK on C' . This 
result corresponds to a test with size a — 0.05 and power 
/3 = 0.44. Requiring that the test have a power /3 of at least 
0.5 relaxes the 95 per cent confidence upper limit to a value 



— 45 ^K. We have tabulated these limits along with 
tnose tor moaels witn coherence scales ot u.O ana z arcmm 
in Table 1. 



7 A COMPARISON WITH OTHER 
MEASUREMENTS OF CMB ANISOTROPY 

The best upper bounds reported to date in observations 
with the Owens Valley radio observatory (OVRO) is that 
the fractional fiuctuations in the CMB sky is limited to 

AT/T < 1.7 X lO"'^ at a resolution of 2 arcmin; /To < 
1.9 X 10~^ in the case of fiuctuations with a GACF and 
a coherence scale of 2.'6 (OVRO NCP experiment: Read- 
head et al. 1989). Sensitive imaging observations with the 
Very Large Array (VLA) at 8.4 GHz have been used to 
place a limit of 2.0 x 10-^ on the fractional temperature 
fiuctuations (AT/T) in the CMB at a resolution of 1 ar- 
cmin (VLA experiment: Partridge et al. 1997). Assuming 
that the CMB anisotropy has a GACF form, an upper limit 
of AT/T ^ 2.1 X 10"^ has been set for a coherence scale 
of ifl based on observations at 142 GHz using a 6-element 
bolometer array (SuZIE experiment: Church et al. 1997). 
The SuZIE observations have also been used to derive a 2-a 
upper limit of Q fiat < 26 at an effective I ~ 2340 (Ganga 
et al. 1997). Among these three observations, the first two 
have sensitivities per resolution element that are compara- 
ble to our ATCA observations. The SuZIE observations have 
covered a larger sky region. The ATCA observations pre- 
sented in this work limit Qfiat to 23.6 /iK at l^ff = 4700 and 
limit C^/To to 1.6 x 10^5 foj. GACF form CMB anisotropy 
with 5c = 1 arcmin. These different experiments attempting 
to measure the CMB anisotropy on arcmin scales have com- 
parable upper limits on random phase CMB anisotropy at 
multipoles / > 1000. 
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sky autocorrelation function C{9) by 



I of CMB anisotropy experiments that are done 
e-dish telescopes with a specified beam switch- 
3, the telescope beam is usually described by a 
on the sky that is normalized to unit volume, 
temperature measured at any sky position is 
mean sky brightness temperature weighted by the 
srn. The antenna temperatures obtained at a set 
itions (that are defined by the beam switching 
■e combined linearly to form an estimate of the 
The 'effective' beam on the sky is the same 
of the 2D telescope beam patterns. The 
ponding to the estimate of the CMB anisotropy 
with the 'effective' beam is the spherical har- 
ion of the 'effective' beam pattern, 
ervations for CMB anisotropy made with Fourier 
telescopes, the sky temperature anisotropy is 
the array through the element primary beam pat- 
temperature anisotropy, attenuated by the pri- 
pattern of the individual array element, is con- 
the synthesized beam pattern. The synthesized 
nc rmalized to peak unity and has zero volume. Un- 
of the 'effective' beam in beam-switched, single- 
the synthesized beam is not usually de- 
e into a linear combination of identical (position 
ams that have finite volume. Therefore, the mea- 
(image pixel intensities) in Fourier synthesis im- 
itegrals of the sky brightness temperature over 
zed beam and cannot be converted to mean (or 
r lean) temperatures by normalising with any beam 
measurements will be in units of flux density 
5am rather than temperature, 
ge multipole orders and small angles, the spheri- 
decomposition may be approximated by a con- 
Folirier transform and the Ci 's may be related to the 



obser /ations. 



Tie 



Ci = 27r / edec{e)Jo{ie). 



(Al) 



Multipole order l/{2n) and angular distance 6 (in radians) 
are Fourier transform conjugates. If the 'effective' beam is a 
purely radial function b{6), the TFF will be similarly given 
by 



F, = 



2n / edeh{e)Jo{ie) 



(A2) 



If the 'effective' beam is not simply a radial function, one 
computes the 2D Fourier transform ]F(Z,^) of the 2D beam 
b(9, (j)) and circumferentially averages |IF(;,0|^ to get the 
fllter function Fi. 

The TFF in Fourier synthesis images will vary across 
the image. At any image pixel location {9o, (j>o), we first de- 
termine the 'effective' beam b{9, (j)) as the product of the 
primary beam centred at the antenna pointing centre and 
the synthesized beam centred at {6o,4>o)- The two beams 
are separately normalized to have peak unity. The beam is 
Fourier transformed to give 



IF(/,C) 



6(0, 0)e 



i{6lcos€^-\-c}>lsin$^) 



The TFF is then derived as 



(A3) 



(A4) 



Fi , as defined in this form for Fourier-synthesis imaging, has 
units (Jy beam~^ K~^)^. 

If the image variance is computed as a linear combina- 
tion of the squares of the pixel intensities (a weighted sum 
of squares) , the TFF for the variance estimator may corre- 
spondingly be computed as the same linear combination of 
the Fi^s evaluated at the individual pixel locations. 
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